How to bosonize fermions with non-linear energy dispersion 
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We develop a systematic method to treat the effect of 
non-linearity in the energy dispersion on the usual bosoniza- 
tion result for the single-particle Green's function of fermions 
in arbitrary dimension. The leading corrections due to the 
quadratic term in the energy dispersion are explicitly calcu- 
lated. In the Chern-Simons theory for half-filled quantum 
Hall systems curvature is shown to be essential. 

PACS numbers: 05.30. Fk, 11.15-q, 67.20+k, 71.27+a 

For many years the bosonization technique has been 
successfully used to study one-dimensional Fermi systems 
beyond perturbation theory [Q. Motivated by experi- 
mental evidence for non-Fermi liquid behavior in high- 
temperature superconductors and quantum Hall systems, 
the generalization of this approach to arbitrary dimen- 
sions d has recently received a lot of attention |^| |J. 
However, although bosonization is non-perturbative in 
the sense that infinitely many Feynman diagrams are 
summed via an underlying Ward identity Jig] , it hinges 
on one essential approximation: the linearization of the 
non-interacting energy dispersion ek- Namely, measur- 
ing wave- vectors with respect to coordinate systems cen- 
tered at points k a (see Fig.pl), we may write £k°+ q = 

2 

+v" q+ 2^ i where v Q is the Fermi velocity and m a 
is the effective mass close to k". Conventional bosoniza- 
tion sets -\ = 0. Haldane [^l) has speculated that 
it should be possible to develop a perturbation theory 
around the non-perturbative bosonization solution of the 
linearized theory, using ^_ as small parameter. An at- 
tempt to construct such an expansion has recently been 
made by Khveshchenko . However, so far his method 
has not been proven to be useful in practice. To the 
best of our knowledge, an explicit calculation of the ef- 
fect of finite m a on the usual bosonization solution for 
the single-particle Green's function does not exist. In the 
present work we shall solve this problem in arbitrary d by 
means of our functional bosonization approach [0 |9) . We 
then use our method to study the relevance of curvature 
in the Chern-Simons theory of the half-filled Landau level 
p3[ . This problem has recently been examined by several 
authors with conflicting results J|,|lJ],[l4|] . We hope that 
our work will clarify the issue. 

For simplicity let us first consider a system of spin- 
less fermions with Landau interaction parameters / kk , 



where q is the momentum transfer between two particles 
with initial momenta k and k'. We start by partition- 
ing momentum space into a finite number of sectors K a , 
a = 1, n, which depend on cutoffs A" and A" as 
shown in Fig.jlj. Assuming that the / kk are dominated 
by |q| ~ q c , we choose A", A" > g c . On the other hand, 
the sectors must be sufficiently small so that the local 
curvature of the Fermi surface is constant within a given 
sector. Furthermore, / kk should not change appreciably 
when k and k' are restricted to given sectors, so that we 
may replace / kk — > f^ a . In particular, for Fermi sur- 
faces with constant curvature and interactions / q that 
are independent of k and k', there is not need any more 
for introducing several sectors as long as we can handle 
the curvature problem. In this case we formally identify 
the entire momentum space with a single sector. 




FIG. 1. Definition of the sectors K a and cutoffs A" and 
\ a . The thick solid line is the Fermi surface, and the dots are 
the origins k a of local coordinate systems. For sectors that 
intersect the Fermi surface we choose k a such that £k» = /i, 
but in general £k° 7^ fJ>- By construction the union [J a K a 
covers all degrees of freedom in the system. 

We are interested in the Matsubara Green's function 
G(fc), which can be represented as a Grassmannian func- 
tional integral in the usual way (lj]|l(J . We eliminate the 
Grassmann fields by means of a Hubbard-Stratonovich 
transformation involving bosonic auxiliary fields <f>" as- 
sociated with the sectors. After the standard transforma- 
tions the Green's function can be exactly written as 
G(k) =< [G]kk >s ef fj where G is an infinite matrix in 
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momentum- and frequency space, with matrix elements where £ Q (q) = e a (q.) — e "(0)> an d u Q (r, r) = 
given by = [G^ 1 ] 



p r $° 



kk> ~ [V]kk>- Here [G ) kk , = 
Skk'Go(k), with Go(fc) = [io; n — £k + /z] — 1 , and the in- 
finite matrix V is defined by [V] kk ' = T, a Q "i k ) V k-k' > 
with Vq = jj<j>q. Here /3 is the inverse temperature, fx is 
the chemical potential, and the cutoff function Q (k) is 
unity if k £ K a , and vanishes otherwise. Let us empha- 
size again that the above construction includes the case 
that we identify the entire momentum space with a single 
sector: then the a-sums run over a single term a — 1 with 
O a (k) = 1, and it is convenient to choose k a such that 
£k= = A*- The symbol < ... >s ef , denotes functional 
averaging with respect to the effective action 5 e // = 

S 2 + S H n, where S 2 = ^ £ 8 E aa >[f^ X } aa ' <t>- g <t>f , and 

Skin = ~Tr ln[l — GoV] . Here V is the volume of the sys- 
tem and f is a matrix in the sector indices, with elements 
-q ' 

[fj aa ' = /q Q - For 9c < A Q , A Q the matrix G is approxi- 
mately block-diagonal, with blocks labelled by the sector 
index a. Shifting k = k" +q and choosing |q| <C A Q , A Q , 
we have 16 1 G(k Q +q, iuj n ) =< [G Q ]qg] >s eff , where the 
Fourier transform Q a (r, r', r, r') of \G a ]qqi satisfies 

[-d r - e Q (P r ) V a (r, t)} G a (r, r', r, r') 

= S(r-r')S*(T-T') . (1) 

Here P r = — iV r is the momentum operator (we use 
units where Ti — 1), e"(q) = £k°+ q , and S*(t) 



jj^2 n e lu>nT is the antiperiodic ^-function. The potential 
V a (r,T) = J2 q e l [ q ' r_a,mT l V q a is the Fourier transform of 
the Hubbard-Stratonovich field V q a = |0« Eq.(|l|) to- 
gether with the boundary condition that a (r, r', r, r') 
should be antiperiodic in r and t' uniquely determines 
CP. Hence, to calculate the Green's function of the 
many-body system, we first need to solve Eq. (|l|) for fixed 
V a (r, t), and then average the result with respect to the 
effective action S e ff. Note that for the action 

S e ff is not Gaussian. The leading non-Gaussian correc- 
tions can be easily calculated dis- 
solving Eq.(Q) for finite m a is more difficult. Let us 
make the ansatz 

g a (r, r', r, t') = Q?{t, r', r, T >) e * a (^)-* a (*' . (2 ) 

To satisfy the boundary conditions, we require that 
<j> a (r, r) should be periodic in r, while (r, r', r, r') 
should be antiperiodic in r and t'. The crucial obser- 
vation is now that we obtain an exact solution of Eq. (Q) 
by choosing <I> Q and Qf such that 



-d T - T(Pr)] * Q (r, t) = V a {r, r) 



[P r $ Q (r, r)f 

2m a 



(3) 



[-9 T -e Q (P r )+/i 



u Q (r,r)-P r ] ^(r,r',r,r') 

= 5(r-r')5*(T-r') , (4) 



Differential equations of the type (3|) are called eikonal 
equations, and appear in many fields of physics, such as 
geometrical optics, quantum mechanical scattering the- 
ory, and relativistic quantum field theories fl7|| . In the 
limit — > the eikonal equation (|^) is linear and can be 
solved exactly via Fourier transformation. The solution 
has first been discussed by Schwinger p8| , see also Refs. 

Furthermore, in this case the velocity u Q (r, r) van- 
ishes, so that Gi(r, r', r, r') = Gq (r — r', r — r'), where 
the Fourier transform of Gg(r,r) is given by G^iq) = 
Go(k a + q, iu) n )- For finite m a , Eq.Q describes the mo- 
tion of a fermion under the influence of a space- and time- 
dependent random velocity u Q (r,r). At the first sight it 
seems that this problem is just as difficult to solve as the 
original Eq.(^). However, perturbation theory in terms 
of the derivative potential u a (r,r) ■ P r in Eq.(|J) is less 
infrared singular than perturbation theory in terms of the 
original random potential V a (r,r) in Eq.([j]). Moreover, 
for large effective mass m a the random velocity u Q (r,r) 
is a small parameter which formally justifies the pertur- 
bative treatment of the derivative potential. 

Let us first consider the eikonal equation (||) . Although 
it is impossible to solve this non-linear partial differential 
equation exactly, we can obtain the solution as series in 
powers of V a jl7| . At this point it is convenient to work 
in Fourier space, 
and = [ioj m 
as 



i(q-r- 



Defining $ a (r, r) = J2 q e 
- £, a {q)]$ni the eikonal can be written 



Q 



<J> Q (r,r) - $ Q (r',T') 



J"(ry,T,T') = G^(q) 



(5) 



e i(q-r-w m T) _ e i(q-r'-w m -r' 



(G) 



where the functional satisfies for q ^ the non-linear 
integral equation 



g q ~ / t<i>i +i ig'q" * q' ^ q" 

q'q" 



(7) 



with the kernel given by ^, q „ = ^rG^(q')G^(q"). 
Here Gg(q) = [ioj m — ^(q)] -1 i s a bosonic Matsub- 
ara Green's function with energy dispersion given by 
the excitation energy £ Q (q) Jl^]. The q = 0-term re- 
quires a special treatment. From the definition = 
[zw m -£ Q (q)]<I>£ it is clear that ^ = 0. Iterating Eq.(Q), 
we obtain a series in powers of the random potential, 
n = J2n=i K, q , where for ^ 



*n, 9 = E S q , qi+ ... +qn C^qi---qn)V q a 1 ---Vql 



(8) 



qi...q„ 



with G" tx (l/m a ) n . The first two vertices are 
Gf ( qi ) = 1 and C%{ qi q 2 ) = 7? V Having solved Eq.(§) 
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to a certain order in V a , we know also the random veloc- 
ity u Q (r, t) in Eq.(^) to the same order in V a . In Fourier 
space Eq.(||) is equivalent with the Dyson equation 



Q.Q' 



■G?{q\q') 



3" 



(9) 



where the matrix elements of the derivative potential are 
®%qi — ^g-q'^,^' wrtn dimensionless vertex \~ ~, — 

Go(q — q') ^ q ^J' q ■ Iteration of Eq.(|9|) generates an ex- 
pansion of §i in powers of the derivative potential. We 
would like to emphasize that we are not simply expand- 
ing in powers of Because the Gaussian propagator 
of the y Q -field is proportional to the screened interac- 
tion ,f PPA ' a within random-phase approximation (RPA) 
J^-^|, the effective expansion parameter is proportional 
to f PPA ' a /m a . This will become obvious shortly. 

To obtain the Green's function of the many-body sys- 
tem, we need to average Eq.(||) with respect to the ef- 
fective action S* e //. Because averaging restores trans- 
lational invariance, we may set r' = r' = and cal- 
culate G Q (r,-r) = (C/ Q (r,0, r, 0)) s . We parametrize 
the average Green's function as G Q (r,-r) = [G"(r, t) + 
G£(r,T)]e« Q ^\ where G?(r,r) = {Qf (r, 0, r, 0)) s 



r(r,r)=ln 



-*<*(0,0)\ 



(10) 



(Sg?(r, 0, t, 0)<Je*°( r ' T )-*°(°'°)) „ 

^r.Tj- / e *«( r ,T)-*«(0,0)\ • 

Here 5X = A— < X >s c}r We now perform the av- 
eraging perturbatively. Note that the fermionic degrees 
of freedom have been completely eliminated, so that the 
perturbation theory is formulated in terms of the bosonic 
field V*. Following Refs. |,|, we calculate Q Q (r,r) via 
a linked cluster expansion. In this way we obtain an ex- 
pansion Q a (r,r) = Yln=l Qn( r > r ): with 



£(r,r) = S q , qi+ ... +qn WZ( qi ...q n ) 
qqi---q n 

x J" (r, t) J"" (r, t) • • • J" (r, t ) 



(12) 



where J qi {p, r) = Jg"(r, 0, r, 0), and the vertices W" can 
be calculated perturbatively in powers of our small pa- 
rameter f RPA ' a /m a . At the leading tree-level (where 
bosonic loops are neglected) the vertex W" is propor- 
tional to (f RPA - a ) n /(m a ) n ~ 1 . Already the first term 
Qi (r, r) contains non-trivial effects due to the finiteness 
of —s. Hence, to study the relevance of curvature, it is 
sufficient to calculate Qf(r, r) at tree- level, in which case 
Wf is approximated by —^vfq 1 PA ' a - This amounts to 
averaging in Eq. _|Hj) with respect to S e f t in Gaussian ap- 
proximation J?]-|9j], and yields Q"(r,r) = R% — 5"(r, r), 
with R<? = S*f(0, 0) and 



q 



/f PA '« cos(q • r - uj m r) 



e(q)]K + e(-q)] 



(13) 



Note that £ Q (-q) = -£ a (q) + ^s, so that for finite m a 
the integrand in Eq.(|13|) has only simple poles. In con- 
trast, for jjjjj = the denominator in Eq.(|l3|) gives rise 
to a double pole, which leads to rather peculiar features 
in the analytic structure of the Green's function in d > 1 
P, p0[. It is also not difficult to calculate the next order 
mf RPA ' a /m a 11. Then one should retain the vertex 
1^(9192) a t tree- level (in which case it is approximated 
by lpbylq 1 q 2 j'q H i PA ' a 'fqY^) and includc one-loop cor- 
rections to W"(qi) (which lead to a small renormaliza- 
tion of the RPA- interaction in Eq.(p"3|)). 

Next, consider the calculation of G" and G%- Because 
the derivative potential removes possible infrared singu- 
larities, we may use the conventional impurity diagram 
technique pi|| . For the average G" —< > we cal- 
culate the irreducible self-energy in self-consistent Born 
approximation. By expanding the self-energy (and not 
directly the Green's function) we take into account an 
infinite number of terms in the iteration of the Dyson 
equation (^|). A truncation at a finite order would lead to 
unphysical multiple poles in the Fourier transform G"(g) 
of G?(r,r). We obtain [Gf(q)}- 1 = [G^(q)}- 1 - E?(g), 
where, to leading order in f PPA ' a /m a , 

S ?® - ^Y,f* PA ' aG o(«')GU-q')^f x 
\^G'l(q + q') + ^[GUq + q')-GUq-q , )] \ ■ 



(14) 



Similarly, we use perturbation theory to calculate 
G%. The leading contribution to Eq.(|ll"|) is of or- 
der f RPA ' a /m a , and yields for the Fourier transform 
Ga (9) = G?(q)Y a (q), with 



Y a (q) 



^G«(q + q') + ^- [GUq + q') - G?(q - q')] • 

I m a m a I 

(15) 

Eqs.(p^|), (|l4| ) and ( |l5| ) are the main result of this work. 
Higher order corrections involve at least an additional 
power of f RPA ' a /m a . For ^ = we have G£ = 
and G" = Gq. In d = 1 we reproduce then the well- 
known bosonization solution for the Tomonaga-Luttinger 
model. Furthermore, direct expansion of our result for 
G Q (r,r) to first order in f PPA ' a exactly reproduces the 
so-called GW-approximation for the self-energy , with 
non-linear energy dispersion. If we had set S" = Y a = 0, 
we would have obtained a discrepancy with the GW self- 
energy, because for finite m a the exponentiation of 
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the perturbation series is not quite correct. In a sense, we 
have exponentiated "too much", so that it is necessary 
to introduce correction terms in the prefactor. However, 
even in d = 1 these corrections can be calculated pertur- 
batively, because the extra powers of q' in the numerator 
of Eqs.(0) and © remove the infrared divergencies. 

Our functional bosonization approach can be general- 
ized to include transverse gauge fields HJ|^], so that we 
have now a powerful non-perturbative method for study- 
ing the relevance of curvature in the Chern-Simons the- 
ory for the half-filled Landau level jD| . To leading order 
in the relevant small parameter (see below), we simply 
need to replace f* PA > a in Eqs.@, @ and @ by the 
corresponding propagator f^ S ' a of the transverse gauge 
field fl^Jf, In the most important parameter regime 
\oJ m \ ^ u jfI*1| ^ v f k (here k — e 2 m*/e is the Thomas- 
Fermi screening wave- vector in d — 2, v* F is the effective 
Fermi velocity, to* is the effective mass, and e is the di- 
electric constant) f^ S ' a can be written as 



/, 



CS,a 



2tt 



1 - (k Q ■ q) 2 



to* \u m \/{v* F \q\) + \q\/q e 



(16) 



where we have used the Coulomb gauge and assumed 
a circular Fermi surface. Here k Q = J^t, q = , 

and Qc = (4>m*vVj^/K, with = 2 |jl3). Substituting 
Eq. ( fig ) into Eq. ( |l3| ) and rescaling the variables J^3| , we 
find Qi(t,t) = gF(f\\ ) rx,f;g), where g = q c /(<j)m*v F ) 
is a dimensionless parameter, fy [f±] is the component 
of the dimensionless vector q c r parallel [perpendicular] 
to k Q , f — v* F q c T, and F(ru,f±,f;g) is a dimension- 
less function that can be explicitly written down as a 
three-dimensional integral. Linearization of the energy 
dispersion corresponds to setting g = in the integrand 
before doing the integration. Then it is easy to show that 
F(f || , 0, 0; 0) ~ In |f|| | for |f|| | -> oo, implying an al- 
gebraic singularity in the momentum distribution [^j , just 
like in a Luttinger liquid. On the other hand, evaluation 
of Eq. ( |l3|) with non-linear energy dispersion corresponds 
to keeping g finite in the integrand. For small g we find in 
this case F(f || , 0, 0; g) ~ —eg for |f m | — * oo, where c > is 
a numerical constant. Thus, the algebraic singularity in 
the momentum distribution is an artifact of the lineariza- 
tion. This gives support to the arguments put forward 
in Ref. pfl , and seems to agree with the experimental 
fact that half-filled quantum Hall systems have a sharp 
Fermi surface [^J. We have convinced ourselves that in 
the one-dimensional Tomonaga-Luttinger model a finite 
value of ^L. does not destroy the algebraic singularity 
in the momentum distribution. Thus, the relevance of 
curvature is a specific property of gauge fields. 

The precise form of the Green's function predicted by 
Eqs.(|f), @ and (jl|) for the half-filled Landau level 
will be discussed elsewhere. Because by construction the 
expansion of these expressions to first order in the in- 
teraction exactly reproduces lowest order perturbation 
theory, the Green's function will certainly not be of the 



Fermi liquid type |l|. For small g it can be justified 
to neglect the higher-order terms in the eikonal expan- 
sion and in the perturbative calculation of E" and Y a . 
Using the estimates of Ref. [Q, we find g ~ 0.6 in the 
experimentally relevant regime |24]. We are thus forced 
to conclude that for an accurate quantitative comparison 
with experiments higher-order terms have to be retained. 
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